Non-local effects in the plasmons of strongly interacting nanoparticles, dimers, and 

waveguides 
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Non-local effects in the optical response of noble metals are shown to produce significant blueshift 
and near-field quenching of plasmons in nanoparticle dimers, nanoshells, and thin metal waveguides. 
Compared with a local description relying on the use of frequency-dependent dielectric functions, 
we predict resonance shifts as large as 10% and field-intensity reduction of an order of magnitude 
at inter-particle distances or metal thicknesses below 2 A. Our results are based upon the specular- 
reflection model combined with a suitable non-local extension of measured local dielectric functions. 
We present a roadmap to design plasmon resonances in nanometer metallic elements with application 
to optical antennas and improved photovoltaic, light-emitting, and sensing devices. 

PACS numbers: 



Knowledge of the optical response of materials at the 
nanoscale has become a pillar of nanophotonics. While 
excellent agreement between theory and experiment is 
generally achieved for nanoparticles and nanostructured 
materials by relying on local, frequency-dependent dielec- 
tric functions [l|, [2J , non-local effects are known to play 
an important role at small distances in the few- or sub- 
nanometer region Q . This range of distances is becoming 
experimentally feasible in metallic dimers nanoshells 
and tips [f|, which exhibit large near-field enhance- 
ment 0, Si with practical application to improved SERS 
biosensing Q , photovoltaics , and solid-state light- 
ning However, no account of non-local effects has 
been reported in such systems, due in part to the com- 
plexity of non-local descriptions of the optical response 
in inhomogeneous environments. 

First-principles calculations of the optical response are 
only available for relatively simple systems, such as bulk 
metals and small clusters of noble metal atoms [12j . How- 
ever, several approximate prescriptions have been elabo- 
rated to deal with more complicated systems, and most 
notably the d— function formalism of Feibelma n 1 1311 and 
the specular reflection model (SRM) 0, [H, lid]. The 
latter is particularly advantageous because it permits ex- 
pressing the response of bounded homogenous media in 
terms of their non-local bulk dielectric function, as al- 
ready reported in studies of fast electrons and ions mov- 
ing near planar surfaces 
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In this Letter we carry out SRM calculations of dimers, 
shells, and thin waveguides made of gold and silver. We 
show that non-local effects produce significant plasmon 
blueshifts and reduction of the near- field enhancement. 
The response is constructed from measured optical con- 
stants, which are extended to include non-local effects in 
the contribution of valence electrons. The validity of this 
model is corroborated by successfully explaining exper- 
imentally observed quantum-confinement effects leading 
to broadening and blueshift of plasmons in small metal 
particles. We draw conclusions on more complicated sys- 



tems of relevance in plasmonics by applying the same 
model. 

A first example of the blueshift produced by non- 
locality is shown in Fig. Q] for scattering spectra of 
gold dimers calculated with (solid curves) and without 
(dashed curves) inclusion of non-local effects. The shift 
becomes more apparent as the inter-particle distance is 
reduced and it is accompanied by sizeable plasmon broad- 
ening. This lowest-order longitudinal dipole mode of the 
dimer involves huge piling up of induced charge in the 
particle junction region |8[, thus enhancing the role of 
non-locality. 

In the SRM used throughout this work, each interface 
separating two homogeneous media is described by sur- 
face charge distributions on either side. These surface 
charges are used to construct boundary fields as if they 
were surrounded by an infinite homogeneous medium. 
The surface charges are determined by imposing the con- 
tinuity of the electric potential and the normal displace- 
ment across the interface 17[. The response of the in- 
homogeneous system is thus reduced to the knowledge 
of non-local bulk dielectric functions of the materials in- 
volved. 

Homogeneous media are conveniently described in mo- 
mentum space q, in which non-locality shows up as a 
dependence of the dielectric function e(q,uj) on q. In 
general insulators behave as local media, whereas delo- 
calized s valence electrons in noble metals constitute a 
source of non-locality. It is thus reasonable to treat the 
response of both insulators and core-polarization in met- 
als as local and to approximate the valence-electron con- 
tribution in the q — * limit by a Drude formula, e D (w) = 
1 — u>p/u)(w + ir)), where uj p is the valence-electron-gas 
plasma frequency and rj is the plasmon width. Delocal- 
ized valence excitations produce a significant q depen- 
dence, as studied by Lindhard in the so-called random- 
phase approximation (RPA) and later by Mermin 
to include electron-motion damping in a self-consistent 
fashion [3]. We adopt here the dielectric function pro- 
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posed by Mermin [e M (q, u>)] to describe the full non-local 
contribution of valence electrons. This function has the 
desirable q — > limit e M (0,uj) = e D (ui). The resulting 
non-local permittivity of the metal is thus approximated 

by M 
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{u)-e u {u) + e M {q,u), 



(1) 



where e loc (w) is the local, frequency-dependent part of 



the response, taken from optical measurements 21] . We 
correct the latter in Eq. {1} by subtracting the local 
valence-electrons contribution e D and adding the non- 
local valence-electrons response e M . Finally uj p and r\ 
are chosen to fit both the long-wavelength tail and the 
observed plasmon frequency of e loc [13] (Jko p = 9 eV and 
hr] = 0.05 eV for gold; huj p = 9.1 cV and hrj = 0.02 eV 
for silver). 

The momentum and distance dependence of the re- 
sponse of bulk gold is represented in Fig. [2j In par- 
ticular, Fig. [2a) shows the so-called loss function of 
gold [Im{ — l/e(q, tj)}] under the approximation of Eq. 
(JXJ) , which is directly accessible to electron energy- loss 
spectroscopy [23}. Non-local effects become dramatic 
in the optical region for q > lnm -1 (i.e., to the left 
of the Landau cutoff for excitation of valence electron- 
hole pairs, u> — qvF, where vp is the Fermi velocity), 
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FIG. 1: Extinction spectra of dimers formed by spherical gold 
particles of 20-nm in diameter. Local (broken curves) and 
non-local (solid curves) calculations are compared for sev- 
eral separations between the particle surfaces, d. Consecu- 
tive curves are offset vertically by 500 nm 2 for clarity. The 
external electric field is parallel to the inter-particle axis. 



volved in plasmonic structures. Actually, the real-space 
response y non - loc obtained from the Fourier transform of 
4ir/[q 2 e(q, u>)] differs considerably from the local approxi- 
mation V loc — l/[re loc (w)] up to distances above ~ 1 nm 
[Fig. dfb)]. The ratio between these two quantities has 
a real part significantly smaller than 1 over that range 
of distances, while the imaginary part lies below 0.2 in 
modulus within the area explored in Fig. [2jb). 
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FIG. 2: (a) Non-local wavelength- and momentum-dependent 
response function of homogeneous gold, (b) Distance depen- 
dence of the non-local interaction potential derived from (a) 
and normalized to the local potential l/[re loc (a>)]. 
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FIG. 3: Extinction spectra of spherical gold particles embed- 
ded in glass (e = 2.3). The present non-local theory (dashed 
curves) is compared to experiment (solid curves, taken from 
Ref. [3|) for various values of the particle diameter (see text 
insets). The calculated spectra have been shifted 15 nm to 
the red to improve comparison. 

It is useful to apply this formalism to spherical par- 
ticles, the optical response of which is conveniently de- 
scribed by their multipolar response to external poten- 
tials of radial dependence r , which develop an induced 
potential — ai/r l+1 , where I is the orbital momentum 
number and u\ is the polarizability of order I. Apply- 
ing the SRM, one finds ^ 
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a is the particle radius, eh is the permittivity of the lo- 
cal host medium surrounding the particle, and ji is the 
spherical Bessel function of order I. 

Fig. [3] shows good agreement because the present non- 
local calculations (broken curves) and experimental data 
Q (solid curves) for light scattering by small gold parti- 
cles. Quite different from the size-independent spectral 
profiles predicted by a local description, non-local effects 
give rise to significant broadening for particle radius be- 
low 2 nm [cf. Fig. [2jb)]. These are the so-called quan- 
tum confinement effects, which show up when the mean 
free path of valence electrons is comparable to the par- 
ticle diameter Q. This agreement between experiment 
and theory supports our approach to describe non-local 
metal behavior at short distances. 

For particle dimers [Fig. [T] we describe each sphere 
by the polarizabilities of Eq. ^ and then perform mul- 
tiple scattering to obtain the self-consistent field 24j, 
including all multipoles with I < 20 to achieve conver- 
gence. The dependence of the lowest-order dipole-mode 
wavelength on inter-particle separation in gold and silver 
dimers is shown in Fig. [4] with (solid curves) and without 
(dashed curves) inclusion of non-local effects. Besides 
significant blueshifts, non- locality increases the plasmon 
broadening (shaded areas), an effect that we attribute 
to the availability of additional loss channels (valence 
electron-hole pairs) in the metallic response. 

Accompanying the far-field phenomenology discussed 
above, the near field is also influenced by non-local ef- 
fects, which can reduce the field enhancement at the 
center of a dimer by nearly an order of magnitude, as 
shown in Fig. [5] This limits the applicability of metal 
structures to SERS and non-linear optics, as compared 
with the expectations opened by the field enhancement 
predicted in the local approximation. Again, losses origi- 
nating in electron-hole-pair excitations contribute to this 
reduction, and so does the finite penetration of the in- 
duced charge towards the interior of the metal by a dis- 
tance ~ vf/oj [13] (see Fig.[5]). This is in contrast to the 
vanishing penetration associated to local descriptions. 

Finally, non-local effects are not exclusive of small par- 
ticles or inter-particle gaps. They are also important in 
other strongly-interacting metallic geometries typically 
encountered in plasmonics like nanoshells [Fig. EH a)] and 
waveguides [Fig. [Sfb)]. Similar to dimers, the plasmon 
modes in these structures undergo bluedshifts exceed- 
ing 10% at thicknesses below 5 A. Interestingly, the mu- 
tually complementary geometries of the dielectric gap 
and the metal film considered in Fig. G^b) should yield 
the same dispersion relation in the non-retarded limit 
(dashed curve), but non-local effects produce different 
dispersion relations for each of these systems. 
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FIG. 4: Dimer dipole-plasmon wavelength as a function of 
inter-particle surface spacing d for gold (a) and silver (b) 20- 
nm nanoparticles. Solid and broken curves correspond to non- 
local and local descriptions of the metal, respectively. The 
width of the modes is represented by shaded regions cover- 
ing the FWHM span of the scattering-cross-section plasmon 
peaks. 



More elaborate methods have also been used to de- 
scribe non-locality, particularly using density-functional 
theory [13] and plasmon quantum models [25], although 
they become impractical in large systems or non trivial 
geometries. In contrast, the SRM used here accounts for 
the main features of non-locality and is sufficiently ver- 
satile to deal with relatively complex systems like those 
discussed above. 

In summary, we have shown that non-local effects are 
relevant in the response of metal systems involving short 
distances below a few nanometers. In particular, we have 
reported large blueshifts and reduced near-field inten- 
sity compared to a local description in metallic dimers, 
nanoshells, and thin films. These effects should be impor- 
tant in the design of plasmonic elements on the nanome- 
ter scale, for which the formalism provided here provides 
a suitable and versatile framework. 

I would like to thank Javier Aizpurua, Jeremy Baum- 
berg, Peter Nordlander, Angel Rubio, and Mark Stock- 
man for fruitful discussions. This work was sup- 
ported by the Spanish MEC (NAN2004-08843-C05-05 
and MAT2007-66050) and by the EU-FP6 (NMP4-2006- 
016881 "SPANS"). 
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FIG. 5: Enhancement of the electric field intensity relative to 
the incident field (\E / E eyit \ 2 ) at the center of metallic dimers 
with the same dimensions as in Fig. [4] ( se e point P in the 
inset of that figure), (a) Spectral dependence of the enhance- 
ment for various separations between the surfaces of 20-nm 
spherical gold particles, (b) Maximum field enhancement as 
a function of separation for gold and silver dimers. The re- 
sults from local and non-local descriptions of the response are 
shown by solid and dashed curves, respectively. 
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FIG. 6: (a) Dipolar mode wavelength of gold and silver shells 
surrounding a silica core (e = 2.1) as a function of metal 
thickness. Solid and broken curves correspond to non-local 
and local descriptions of the metal, respectively. The width of 
the modes is represented through shaded regions covering the 
FWHM span of the scattering-cross-section plasmon peaks, 
(b) Plasmon dispersion relations in gold-silica-gold and silica- 
gold-silica waveguides for a 0.5-nm thick central layer (solid 
curves); a non-local description yields the dashed curve in 
both cases. 
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